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Abstract 



The level of a function / on R" encloses a region. The volume of a region between two such levels 
depends on both levels. Fixing one of them the volume becomes a function of the remaining level. 
We show that if the function / is smooth, the volume function is again smooth for regular values of 
/. For critical values of / the volume function is only finitely differentiable. The initial motivation 
| for this study comes from Radiotherapy, where such volume functions are used in an optimization 

process. Thus their differentiability properties become important. 

-I— > 
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1 Introduction 

The volume of a set enclosed by two different level sets of a function / : M™ — * R depends on both 
' levels. Here we fix one level and the question we want to address is the differentiability of the volume 

as a function of the varying level for arbitrary dimension n. It will turn out that under mild conditions, 
among which smoothness of /, this function is again smooth for all regular values of /, but only finitely 
differentiable at critical values of /. Moreover we also consider / on a compact subset V of its domain and 
subsequently pose the same question for level sets restricted to V. Then we get a similar differentiability 
result when we also include the levels of / restricted to the boundary of V. In section 2 we will more 
precisely define the situation we consider and give precise statements of our results. 

The initial motivation for this study comes from Radiotherapy. A patient is treated with ionizing radiation 
causing energy release per unit mass or volume inside the patient, which is called dose. The above 
mentioned function / represents this dose and the set V represents a patient's organ or a tumour region. 
The therapeutic outcome of Radiotherapy treatment not only depends on the dose on the tumour but also 
on the dose on healthy organs. Since it is usually hard to estimate the effect of dose in three dimensons, 
one reduces it for each tumour and organ to a so called dose-volume histogram, see [2] . In it one records 
for each dose value which fraction of the volume receives at least this dose. The resulting collection of 
dose- volume pairs is called the cumulative dose-volume histogram, which is closely related to the volume 
function introduced in the first paragraph. It is well established that the larger the fraction of the 
volume of a tumour receiving a prescribed dose, the larger the probability that the tumour is eradicated. 
For healthy organs and tissues the situation is less clear but some evidence exists that damage can be 
estimated from dose-volume histograms, see [2] and references therein. 

In Radiotherapy treatment planning one tries to find an optimum for a sufficiently high dose on a tumour 
and a sufficiently low dose on healthy organs. The object function of this optimization process depends in 
particular on certain dose-volume pairs (d,v(d)) on the graph of the cumulative dose-volume histogram, 
see [5] . In the optimization process we have a family of dose functions and thus a family of dose- volume 
histograms, parametrized by the optimization variables. Anticipating the result on finite differentiability 
of the dose- volume histogram at a critical value of the dose function, we conclude that the object function 
is finitely differentiable at values of the optimization variables for which d is a critical value of the dose 
function. Because many numerical methods to solve such problems assume differentiability of the object 
function to some order, differentiability of the dose-volume function becomes important. 

The remaining part of this article is organized as follows. In section 2 we give a summary of the results, 
preceded by definitions. In section 3 we sketch the idea of the proof and the results are proved in a series 
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of propositions. Some of the more elaborate computations are summarized in an appendix. 

2 Statement of results 
2.1 Definitions 

Our main object is a smooth positive function / which is bounded and whose levelsets are compact. 
These properties are sufficient for our results, we do not claim necessity. We now define a function class 
for future reference. 

Definition 1 Let C be the class of functions f : R" — > R satisfying the following 

a) f is C°°, 

b) f is positive, 

c) f is decreasing, that is for each e > there is a compact K C K™ such that for all x e R" \ K, 
fix) < e. 

We define the levelsets N c of / e C and sets U c enclosed by levels in a straightforward manner. 

N c = {x e R™ I f(x) = c}, 
U c = {x e R n | f(x) > c}. 

We wish to study the study the differentiability of the volume vol n (U c ) of U c as a function of c. Here 
vol n (A) is the standard volume (Lebesgue measure) in R™ of a measurable set A C R™ 



vol n (A) = 




When we restrict / to a subset V of R™ we make similar definitions. To be more precise, let V C R™ 
be a compact, connected, n-dimensional C°°-manifold with boundary dV which is a (n — l)-dimensional 
C°°-manifold. The restriction of / to V will be denoted by We define V c — U c r\V and the volume 
function we consider for this case is vol (V c ). 

As already indicated in the Introduction, the function vol need not be continuous. Indeed, let / be 
constant on an open neighbourhood of V then it is easily seen that vol is similar, upon scaling and 
translation, to the Heaviside function. In order to avoid this and other degeneracies we assume the 
following non-degeneracy conditions on /, f\y and V. 

Definition 2 Non- degeneracy conditions. 

a) Critical points of f are non-degenerate, that is if x is a critical point of f then det(Hess u{x)) ^ 0. 

b) If x is a critical point of f then x ^ dV . 

c) Critical points of f\v are non-degenerate. 

d) f is fine, that is if x and y are critical points and x ^ y then f(x) ^ f(y)- 

e) f\v is fine. 

Remarks 

1. Conditions a, b and c are essential for our proofs. If one of them is not satisfied our standard 
forms, see section 3, of / and V are no longer valid. In order to proceed we would need higher 
order information about / and V. Moreover the critical point would not be stable under small 
perturbations. At the moment we are not interested in such a situation. Conditions d and e are for 
convenience only and can easily be dropped. 

2. In principle, the sets Uh, V and Vh need not be connected. If one of them consists of several 
components, the construction in the following sections can be carried out for each component 
separately. Therefore without loss of generality we may as well assume connectedness. 
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3. In Radiotherapy one sometimes uses the so called di fferential dose- volume histogram, see [2]. How- 
ever in general this is not a function. In terms of our function / it is in fact the image measure of 
the standard measure on R 3 under / on K, see [1]. Instead of studying the volume function via / 
as defined in this section one could also study the volume function via this measure. The present 
approach however seems to be simpler. 

2.2 Results 

In order to state the results we will make a distinction between regular values and critical values of /|y. 
In case of a critical value of f\v we make a further distinction whether the critical point is in dV or not. 
Results for / and f\y are identical when the critical point is not in dV, therefore they are not stated 
separately. Our proofs are valid only for / in class C, see definition 1, satisfying the nondegeneracy 
conditions in definition 2. 

The first theorem states that for all regular values of the function /, which means for almost all values 
by virtue of the nondegeneracy conditions, the volume function is a smooth function of the level. 

Theorem 1 Let e M be a regular value of f\ v then vol(V/j) is a smooth function of h at 0. 

The second and third theorem state that at a critical value the differentiability of the volume function 
is finite. The order of differentiability depends on the dimension of the domain of /. The nature of the 
discontinuity depends on the Morse index of / at the critical point. Together with the proofs we give 
details about the discontinuity in the next section. 

Theorem 2 Let £ R be a critical value of f\v and let £ dV be the critical point. Then the \?f\-th 
derivative ofvol(Vh) is discontinuous at h = 0, all lower order derivatives are continuous. 

Theorem 3 Let € M. be a critical value of f\dv and let € dV be the critical point. Then the l 21 ^-] -th 
derivative ofvol(Vh) is discontinuous at h = 0, all lower order derivatives are continuous. 

In a Radiotherapy setting, which originally motivated this study, the dimension of the domain of /, 
representing dose, is 3. This means that at a critical point of dose, the volume function is not even twice 
continuously differentiable. As mentioned in the introduction one aimes at finding an optimum for a 
sufficntly high dose on the tumour and a sufficiently low dose on healthy organs. The objective function 
in this optimization process depends on the volume function. Many iterative optimization methods use 
a quasi-Newton method in the background and thus require differentiability to second order. Such a 
method is not guaranteed to be well behaved near a critical value of dose. 

3 Proof of results 

We will work in the class of C°°-functions. This class is closed under the action of the group of C°°- 
transformations. Therefore we have the notion of C°°-equivalence of functions. We use this to put the 
function at hand into a suitable standard form. In general however, this standard form is only valid 
locally in a small open ball. 

Now our aim is to compute the volume of the set enclosed by two levelsets of a function, which is a 
global rather than a local problem. A natural way to look at this set is by "sweeping out" using the 
gradient flow of the function. Using compactness we can then turn this into a local problem considering 
flow boxes of the gradient flow starting in small subsets, here we take (n — I)-simpliccs, of the levelset 
corresponding to the lower level and ending in the levelset corresponding to the higher level. By a suitable 
C°°-transformation we turn each flow box into a Cartesian product of a (n — I)-simplex and an interval. 
This greatly simplifies finding the volume of the enclosed set and its dependence on the level. 

A complication in this procedure is that a general C°°-transformation does not map a pair (/, grad/) 
into a new pair (g,gradg). In order to achieve the latter, the transformation would have to preserve the 
inner product which is used to define grad. For our purposes this is too stringent a restriction. Another 
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reason not to use the gradient flow is the following. We also wish to consider the set enclosed by lcvclsets 
restricted to a set V. The gradient flow is not necessarily tangent to dV, thus the flow box we construct 
might not be restricted to V. Instead we use the flow of a C°°-vcctorficld which is only transversal to the 
levelsets and tangent to dV if necessary. These properties are preserved by a general C°°-transformation. 

The tansformations we apply to a flow box do not in general preserve its volume. Here, however, we 
are only interested in the h dependence of this volume, not in its numerical value. Therefore we may 
apply affinc transformations without any further considerations. With other transformations we have to 
be more careful and we will take them into account at the appropriate places. 

Let us sketch the steps in the proofs of theorems 1, 2 and 3. If is a regular value of /, all points in 
Vh a- r e regular for h small enough. We first construct a finite number of boxes Bi covering Vh using a 
triangulation of Nq and a regular flow from Nq to Nh- Assuming the boxes Bi are small enough we put 
each of them in standard form by several local C°° -transformations. The first transformation parallellizes 
the flow from Nq to Nh- The second transformation is linear and preserves the parallellity of the flow 
but makes it perpendicular to A at and parallel to the last basis vector e n of R". The third transform 
takes / into a local standard form preserving all of the previous. The result is that Bi is transformed to 
the Cartesian product of a simplex in No and the interval [0, h]. The conclusion is that the volume of Bi 
is a smooth function of h. 

If is a critical value of / we use the fact that is the only critical point on Nq. We construct one special 
box Bo containing and away from we use the same construction as above. Differentiability is then 
determined by vol(-Bo)- Again we use a transformation that takes / into standard form, but now at the 
critical point. 

The main part of the proof of theorem 1 is the construction of the boxes and putting them into a standard 
form. In the proofs of the other theorems the emphasis is on computing the volume of box Bo ■ 

Where necessary we assume the existence of a standard basis and a standard inner product. 
3.1 Proof of theorem 1 

Proposition 4 Let be a regular value of f, then for sufficiently small h > there exists a finite 
collection of sets Bi with i G I C N satisfying 

a) V h = UieiBi 

b) vol n (V h ) = ^ ieI vol n (Bi). 

Proof of proposition 4. 

Construction. If is a regular value of / then No is a smooth manifold. The non-degeneracy conditions 
imply that critical values of / are isolated, therefore an h > exists such that all h G [0, ho] are regular 
values. Then all Nh are diffeomorphic to No, sec [3]. First we assume that Ao does not intersect the 
boundary dV of V. Since A is a compact C°°-manifold it allows a finite triangulation with (n — 1)- 
simplices cr^n-i and i G / C N, sec [4]. One of the properties of a triangulation is that for i ^ j either 
Cj,n-i H Pj,n-~i = or Gi t n-\ l~l Cj.n-i = Ofc,n-2 for some k and (n — 2)-simplex <Jk,n-2- Let A be a 
C°°-vector field transversal to Nh for all h G [0, ho). More precisely we impose the condition that there 
is an e > such that \{nh{x), X (x)) | > e for all x e Nh and h G [0, ho], where nh(x) is a unit normal 
to Nh at x. Such a vectorfield exists, for example grad/. If A intersects dV then the non-degeneracy 
conditions imply that the intersection of Nh and dV is transverse for all h G [0, ho] . We restrict to A n V 
which is still a compact C°°-manifold. Now we impose one more condition on the vector field A namely 
that it is tangent to dV n V ho . Let $ G Diff(M") be the flow of A with $(x, 0) = x. Finally we define 

Bi = {$(x, t)\x€ o-i, n - U t G [0, T]} n V h . 

Proof of Vh — UigjJBj. Clearly U i£ iBi C Vh- Suppose x G Vh, then 4rf($(x,t)) is strictly increasing 
or decreasing since \{n h {x) , X (x))\ > e. In either case a finite to exists such that /($(x,i )) = 0. This 
means xo = $(x, to) G A , therefore an i exists such that xo G CTi,n-i which implies x G Bi. Compactness 
of Vh guarantees that a T > exists such that for all x G Vh, |to| G [0, T]. The conclusion is that 
Vh = Ujg/Bj. 

Proof of voln(Vh) = J2i£i V °ln(-Bj). It suffices to show that for i ^ j, vol„(Sj n Bj) = 0. Suppose 
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Bi n Bj and x £ Bi n Bj, then there is a to such that xo = <&(a;,to) £ A*o. By definition <&(a;,t) <E 
Bi H Bj therefore xq E Bi D Bj which means xq £ <Ji. n -i H Oj.n-i = crfc jn _2- From this we conclude 
Bi n Bj = {<J>(x,i) | x e (Tk,n-2,t S [0,T]} n V^. But then dim(Bj nBj) = n - 1 and therefore 
vol„(B 4 n b/) =0. ■ 

The next step is to put the boxes Bi of proposition 4 into a standard form. To do this we also need a 
local standard form of the function /. 

Proposition 5 Let f be a function as in definition 1, satisfying the non- degeneracy conditions in defi- 
nition 2 and /(0) = 0. We distinguish three different cases: 

a) is a regular point of f , 

b) is a critical point of f , 

c) is a critical point of f\dv- 

Then an open neighbourhood O of and a diffeomorphism $ exist such that F = takes one of the 
forms: 

a) F(£,ri) = T], with (£,77) e (R"- 1 xR)nO, 

b) F{£, 7?) = ELi (f ~ EiLi Vl with (£, ry) e (R p x M 9 ) n O and p + q ~ n, 

c) F(£, V,0=C + ELi f< - Ei=i % 2 and dv ls 9 wen by( = with (£, 77, C) G (R p xPxljflO and 
p + q = n — 1. 

Remarks 

1. Case b) of proposition 5 is called the Morse lemma. A critical point in this case has Morse index 
q, but some times it is more convenient say it has Morse type (p, q). 

2. Due to non-degenracy condition 2, is a critical point of f\gy as soon as is a critical point of f\v- 
In the proof of 5 we will need a lemma which we only state here, for a proof see [3]. 

Lemma 6 Let f be a C°° function on a convex neighbourhood O C R™ of with /(0) = 0. Then 
f( x ) = Yh=i Xifi(x) for certain C°° functions with /j(0) = ^f-/(0). 

Proof of proposition 5. 

a) Since is a regular point of / there is a nonzero vector a such that grad/(0) = a. After an 
orthogonal transformation we may assume that with respect to coordinates (x,y) £ R" _1 x R we 
have grad/(0, 0) = (0, \a\). Now we define new coordinates by the diffeomorphism $ : (x,y) 
(x, f(x, y)), then F = <&*/ takes the desired form. 

b) See [3]. 

c) If is a critical point of fv, the the tangent spaces of N and V at coincide. By the non- 
degeneracy conditions is a regular point of /. Using the arguments of case a) we assume that we 
already transformed to coordinates (x,y) £ R n_1 x R such that grad/(0,0) = (0, \a\). Furthermore 
we made the assumption that dV is a smooth manifold, so at least locally it is the level set of a 
smooth function g. Now we apply part a) to bring g into standard form, then on new coordinates 
(u, v) £ MJ 1 ' 1 x R, dV is locally given by v = 0. After scaling in the v direction the function / 
satisfies: /(0,0) = 0, gf-/(0,0) = for i £ {1, . . . , n — 1} and ^j/(0,0) = 1. The remainder of 
the proof is only a slight adaption of the proof in [3] for case b), but included here for the sake of 
completeness. Applying lemma 6 to / and its partial derivatives we get 



f(u,v) = ^2 UiUjaij(u,v) + vf n (u,v) 

i,j<n 
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where and /„ are smooth functions and / n (0, 0) — 1. As a first step we apply the transformation 
(u, v) i ► (u, vf n (u, v)) — (u, z), preserving the standard form of g since / n (0, 0) ^ 0. Then / takes 
the form 

f(u,z)= UiUjOiij(u,z) + z 

i,j<n 

with new functions a.^ . From now on we only apply transformations of the form (u, z) i— » (<fi(u, z), z) 
and we proceed by induction. We assume that 

f(u, z) = ±u\ • • • ± u^-i + Y UiUjOtij (u, z) + z 

i,j=k 

for a certain k > 0. Now let in = Ui for i ^ k and z = z and 

n-1 

Uk = y/a kk (u,z)(u k + Y u i a ik(u,z)/a k k(u,z)) 
i=fc+l 

In order to define u k it may be necessary to permute the rows of onj so that a kk (u, z) ^ 0. Such a 
permutation exists because det(Hess/(0, 0)) ^ which means there is at least one i £ {k, . . . , n— 1} 
such that a,fc(0, 0) ^ 0. Then by continuity there is a neighbourhood of (0, 0) such that atik(u, z) ^ 
0. This means that in each induction step the neighbourhood on which our result holds might 
shrink. Since we only need a finitie number of steps this does not cause any problems. Dropping 
the hats we get in new coordinates 

n-1 

f(u, z) — ±uf ■ ■ ■ ± u\ + UiUjCtij (u, z) + z 

i,j=k+l 

Renaming the variables we arrive at the desired form of F = <&*/, where $ is the composition of 
the transformations in each induction step. ■ 



Proposition 7 Let {Bi} ieI be a collection of boxes as in proposition 4- Let x be a point in the interior 
of o~i,n-i f or some i e L, without loss of generality we assume that x = 0. Then a difjeomorphism <f> 
exists such that 



vol n (B l )= [ ldx= [ [ 

J Bi JO J a 



'CT«,n-l 



dn, 



where J$ is the Jacobian o/$. Moreover vol n (Bi) is a smooth function of h. 



Proof of proposition 7. The vector field in the construction of box Bi has no stationary point therefore 
it can be parallellized by a diffeomorphism <j>i, see [4]. The next two transformations preserve parallellity 
because they are linear. By a linear diffeomorphism $2 we can arrange that X is perpendicular to Nq in 0. 
By another linear diffeomorphism $3 we rotate such that X is parallel to the last basis vector of R™ . The 
last diffeomorphism $4 : (x,y) (x, f(x, y j) takes the function / into local standard form, see proposition 
5. Since $4 is a position dependent shift in the direction of the vectorfield, parallellity is again preserved. 
However the parametrization of the integral curves will change in general. Then $ = $4 o $3 o <£> 2 $1 is 
again a diffeomorphism and on new coordinates (£, 77) we have Bi = <7j jn _i x [0, h}. Since J$ is a smooth 
function and vol n (i?j) depends on h only via the upper limit of the outer integral, vol„(Bj) is a smooth 
function of h. ■ 

Using the previous propositions we are able to prove theorem 1. 

Proof of theorem 1. Construct boxes Bi as in proposition 4. Then by the same proposition vol n (V/j) = 
J2iei v0 ^n(Bi). In the latter vol„(Bj) is a smooth function of h by proposition 7. The finite number of 
boxes guarantees that also vol n (Vh) is a smooth function of h. ■ 
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3.2 Proof of theorem 2 

In the previous section all points were regular points of /. Here too points will be regular except the point 
which now is a critical point. Therefore we will use the same construction of boxes as in the previous 
section. Only the box Bo containing will be treated differently. This means that differentiability in this 
situation is determined by vol„(-Bo)- 

A critical point with Morse index is a minimum of /. If the Morse index is n the critical point is a 
maximum of /. Differentiability for minima and maxima is very similar so we only state a result for one 
of them. A critical point with Morse index q where < q < n is called a saddle. 



Proposition 8 Let be a non-degenerate critical point of f with Morse index 0, so is a local mini- 
mum. Ther 
continuous. 



mum. Then the \%]-th derivative o/vol(V^) is not continuous at h = 0, all lower order derivatives are 



Proof of proposition 8. Let h > 0, then is the minimal value of / so V-u is empty. The level set 
N only contains 0. Using the neighbourhood O of proposition 5 on which / takes its standard form we 
compute vo\ n (Vh) using polar coordinates (r, ip), then that u(r, <p) = r 2 . Here we assume that h is small 
enough so that Vh C O. The Jacobian of the transformation <E> in proposition 5 will be denoted by J$ 
and the Jacobian of changing to polar coordinates by r n ~ 1 g n ((p). Since $ is non-singular and C°° we can 
split J$ = c + K$ where c ^ and both J$ and are C°° . Then we have 

vol„(V4) = / ldx= f ,hd^ = cf f r n - 1 drg n ((p)dip+ [ d£. 
Jv h Jv h J Jo Jv h 

For a non-differentiability result it is enough to consider the first integral in the last expression, since this 
integral contains the lowest order terms in h and we are interested in h — > only. Let a n be the "area" 
of S™- 1 then 

c J r n - x dr g n {<p) dip = ca n r^ 1 dr = ^h n ' 2 . 
Thus we obtain the result that the [fl-th derivative is discontinuous at h = 0. ■ 

Proposition 9 Let be a non-degenerate critical point of f . Assume the Morse type of is (p,q) with 
p^O and q^O. Then the \-^~\-th derivative of vo\{Vh) is discontinuous at h = 0. For both p and q even 
the discontinuity is a jump, for p and q odd it is a log-like singularity and for p + q odd it is a root-like 
singularity. All lower order derivatives are continuous. 

Proof of proposition 9. Let h be small enough so that is the only critical point of / in V4 - 
Furthermore let O be a neighbourhood as in proposition 5 such that / can be put into standard form c) . 
Changing to cylinder coordinates (r, ip, s, ip) we may assume that f(r, ip, s, ip) = r 2 — s 2 . Then an e > 
exists such that B = {(r, tp, s, ip) £ M x 5 P_1 x R x | r + s < e} is a subset of O. The intersection 

of dB and TVo is transversal so again taking ho small enough we may assume that dB transversally 
intersects Nh for all h e [0, ho]- Then Vh \ B is a compact C°°-manifold (with boundary) so we can apply 
the construction of boxes as in proposition 4. Where we have to impose the additional condition that the 
vector field X is tangent to the boundary dB of B. By proposition 7 J2 ieI vol„(Sj) is a smooth function 
of h. If we now set B = B n Vh then vol„(V; l ) = X^e/ V °l™(^) + v °ln(-Bo)- Thus differentiability of 
vol n (Vft) is determined by the differentiability of vol„(_Bo). 

The Jacobian of the transformation $ in proposition 5 will be denoted by J$ and the Jacobian of changing 
to polar coordinates by r p_1 s 9_1 g P ((p)g q {ip) ■ The last transformation we apply is a scaling so that B is 
bounded by r = 0, s = and r + s = 1. Then we get 



vol„(.Bo) = / ldx= J^d^dr/^c ldS,dr/+ K^d^dn. 

<J Bq <J Bq <J Bq <J Bq 



For our result we only need to compute the first integral in the last expression, using the same arguments 
as in the proof of proposition 8. The actual computation can be found in appendix A.l. From lemma 12 
in the same appendix the result on differentiability follows. ■ 

Proof of theorem 2. The proof follows from propositions 8 and 9. ■ 
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3.3 Proof of theorem 3 

In this section is a critical point of f\ov, but all other points are regular. First we consider a critical 
point which is a minimum of f\gv- Since for a maximum we get the same result we do not state it 
seperately. 

Proposition 10 Let be a critical point of f\gy. Assume the Morse type is (n — 1, 0) with n > 1. Then 
the |~^±i]-s< derivative of vo\(Vh) is not continuous at h = 0, all lower order derivatives are continuous. 

Proof of proposition 10. The proof is very similar to that of proposition 8 therefore we will only 
indicate the essential differences. On a neighbourhood O of such that / can be put into standard form 
c) of proposition 5 we take cylinder coordinates (r, ip,Q elx R"~ 2 x R. Then locally f(r, <p,Q = r 2 + ( 
and dV is given by ( — 0. Splitting the Jacobian as in proposition 8 we obtain 



vol n (V h )= ldx= .Udi = c / r n - 2 d(drg n _ 1 (^)d^+ K^d^. 

JV h JVh J JO JO JVk 



rvh rh- 

>V h JV h J JO JO 

Only evaluating the first integral (cf. proof of proposition 8) we get 

r-Vh rh- 



J JO JO 



~"- 2 dCdrg n ^) d V = 2c^-/ i ("+ 1 )/ 2 . 



Thus we obtain the result that the [^y^-st derivative is discontinuous at h = 0. ■ 
Next wc turn our attention to saddle points of f\av- 

Proposition 11 Let be a non-degenerate critical point of f\ov ■ Assume the Morse type of is (p,q) 
with p ^ 0, q ^ and p + q = n — 1. Then the f 2 ^-] -st derivative of vol(Vft) is not continuous at h = 0. 
For both p and q even the discontinuity is a jump, for p and q odd it is a log-like singularity and for p + q 
odd it is a root-like singularity. All lower order derivatives are continuous. 

Proof of proposition 11. We proceed along the lines of the proof of proposition 9 again indicating the 
main differences only. Let O be a neighbourhood of such that / can be transformed to standard form 
c) of proposition 5. Taking cylinder coordinates (r,<p,s,ip, q() we may assume that f(r,ip,s,rh,() = r 2 — 
s 2 + (. Then an e > exists such that B = {(r, tp, s, i>, C) <E Rx S' 9 ^ 1 xlxS'^'xR | r + s < e, < C < e} 
is a subset of O. Once again taking ho small enough we may assume that dB transversally intersects 
for all h £ [0, ho]. Then Vh \B is a compact C°°-manifold (with boundary) so here too we can apply the 
construction of boxes as in proposition 4 with the additional condition that the vector field X is tangent 
to dB. We set Bo = BP\Vh then vo\ n (Vh) — J2i£i v0 ^n(Bi) +vol n (So). Now differentiability of vol„(Vft,) 
is determined by the differentiability of vol„(i?o)- Splitting the Jacobian as in proposition 9 we get 

vol„(B ) = ldx= J$ d£ drj d( = c ld£dnd(+ K<s>d£dn d(. 

J Bq J Bq J Bq J Bq 

We only compute the first integral in the last expression (cf. proof of proposition 8). For the actual 
computation see appendix A. 2. From lemma 13 in the same appendix the result on differentiability 
follows. ■ 

Proof of theorem 3. The proof follows from propositions 10 and 11. ■ 
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A The volume of Bq containing a saddle point 
A.l The volume of B in proposition 9 

We first recall the definition of Bq. Let O be a neighbourhood as in proposition 5 such that / can be 
put into standard form c). Changing to cylinder coordinates (r, <p, s, ip) and after an appropriate scaling 
we may assume that f(r, ip, s, ip) = r 2 — s 2 and B = {(r, ^s^jelx 5 P_1 x M x S* 9-1 | r + s < 1} is 
a subset of O. Then we define B n = {(r, <p, s, ip) e B | < r 2 — s 2 < h} where h > 0. Now we wish to 
compute 



/ ld£dri= / r p 1 s q 1 g p (4>)g q {ip)dr ds d(f>dip = Cp, q I p ,q{h). 

J B •> B 



The Jacobian of the transformation to cylinder coordinates is given by r p_1 s 9_1 g P (4>)g q (ip) ■ The integrand 
does not depend on the angles so we split off the angular part and denote the integrals by a p and a q where 
a m is the 'area' of the m — 1 sphere. To facilitate the computations we make one further transformation: 
u = r + s, v = r — s. In the constant c p _ q we absorb the constants a p , a q and the Jacobian of the change 
of coordinates (r,s) to (u, v). Furthermore we distinguish h < and h > and to simpify notation we 
write k(u, v) = (u + w) p_1 (m — See figure 1 for the regions of integration. 



W ft ) = Iv=hl-u k(u,v)dvdu, h<0 
Ip.q {h) = Jo f-u k( u > v )dv du + j Q u k(u, v)dv du + k(u, v)dv du, h > 




After some computations it turns out that I p , q (h) consists of several parts: for both h < and h>0: 
fi q (h) = P£ q {h) + a± q {±h)^' 2 + f3± q (±h)^ 2 logV±h + 7 ±,. 

There is a constant part because we consider the volume between levels —1 and h. This is done for 
computational reasons. Then continuity demands that ^f p q — "f p q which turns out to be true. Common to 
all cases is a polynomial part P p q (h). The part with coefficient a p q may contain a square root depending 
on p and q and if the coefficient [3 p q is nonzero there is a logarithmic part. The definitions of P, (3, 7 
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and a are as follows 

p+(h) = p-(h) = V () y k ){ m ){ L j h 

v«\ ) p,^ > - {k + m+l){p + q-2{k + m + l)Y 

k.m 

a+ =B - = V ( 2 )(V)( 9 ™)(-ir 



(p + q) 

^PA ~~ ^P,Q ~ 



k,m 



hm (p + q)(p + q-(k + m+l)y 



ct^ = — a 



^ (p + g )(p + g -2(fc + m+l)) ^ + 
2 2 



(p + g) v y b + ?) 2 

The sum Y^kln IS taken over all k and m satisfying < k < p—1, < m < q— 1 and 2(fc + m+ 1) 7^ (p+<?), 
whereas the sum Y^km 1S taken over the same range of k and m but now 2(k + m + 1) = (p + q). The 
last line defines cti and C7 2 . With these definitions (3+ q = ^ q) o- 2 (p, q). The following properties of o\ 
and (72 are easily checked. 

1. (Ji(p,q) = (-l) q ai(p,q), so ai(p,q) = for q odd, 

2. a 2 (p,q) = (-l) q - 1 o 2 (p,q), so a 2 (p,q) = for q even, 

3. (T2(p, <?) 7^ only if p + q even, so using 2) 02 (p, 9) 7^ only if both p and q are odd. 
Now the next lemma is immediate. 

Lemma 12 For each p and q the \^~\ -th derivative of I Piq as a function of h is discontinuous at h = 0. 
The nature of the discontinuity depends on p and q. For both p and q even it is a jump, for p and q odd 
the discontinuity is a log-like singularity and for p + q odd it is a root-like singularity. 



A. 2 The volume of B in proposition 11 

The computation of the volume of Bo in proposition 11 is similar to that in section A.l. Here we only 
indicate the differences. 

First note that p and q have a slightly different meaning because p + q = n — 1. In this case the function 
k in the expression for I p , q (h) is given by k{u, v) = (h — uv){u + v) p (u — With this definition of 

k the functions I^ q (h) are defined as before. Again after some computations we find 

Ip, g (h) = P± q (h) +a p ± g (±/ l ) (p+9+2)/2 +/? p ± g (±/ l ) (p+9+2)/2 logV±^ + 7p ± 9 + 8± q h. 

The expressions for P, a, (5 and 7 are more involved than in the previous section. Their structure, 
however, is similar therefore we skip the details. The relations are equal. There is one new term in the 
expression above which is defined as 

p - q p ' q " (p + q)(p + q-(k + m+l)y 

Thus we come to the same conclusion as in section A.l. 

Lemma 13 For each p and q the \ !] 4r-~\-st derivative of I p , q as a function of h is discontinuous at h = 0. 
The nature of the discontinuity depends on p and q. For both p and q even it is a jump, for p and q odd 
the discontinuity is a log-like singularity and for p + q odd it is a root-like singularity. 
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